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Superfluid "^He-A has domain-waU-like structures, which are called solitons. We calculate nu- 
merically the structure of a splay soliton. We study the effect of solitons on the nuclear-magnetic- 
resonance spectrum by calculating the frequency shifts and the amplitudes of the soliton peaks for 
both longitudinal and transverse oscillations of magnetization. The effect of dissipation caused by 
normal-superfluid conversion and spin diffusion is calculated. The calculations are in good agree- 
ment with experiments, except a problem in the transverse resonance frequency of the splay soliton 
or in magnetic-field dependence of reduced resonance frequencies. 
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' Nuclear magnetic resonance (NMR) has turned out to be very useful for studying the superfluid phases of liquid 
(2jr): '^He. The two superfluid phases A and B are distinguished in the NMR spectrum by different frequency shifts of the 
^ absorption peaks. ^ In addition to these "bulk" peaks, one often observes additional "satellite" peaks. These are caused 
by topological objects and textures that appear in the superfluid order parameter. Especially in superfluid ^He-A, 
several different objects have been identified based on the frequency shifts of the satellite peaksn^iii The simplest of 
these are solitons. They are domain-wall-like structures where a planar object separates two different but degenerate 
bulk states. 

■ The satellite peaks in ^He-A were first observed in measurements in the mid 1970's,-'^ Soon after the theory of 
^ , solitons in '^He was developed by Maki and Kumar Their calculation gave a striking agreement with the measured 

frequency shifts of the satellite peaks at temperatures close to the superfluid transition temperature Tc- This initial 
•/"i , success had the consequence that further studies of solitons went to other directions^iiSiiiiiSi and, unfortunately, no 

■ more precise calculations were done. 

\ A soliton has two basic structures, "twist" and "splay," which correspond to the cases of a magnetic field {B ^ 1 

■ mT) perpendicular and parallel, respectively, to the plane of the soliton wall. Both these structures can be studied 
' using small oscillations of the magnetization that are either transverse or longitudinal relative to the static field. 

"Y"; There are several points that can be improved in the previous calculations, given as follows, (i) The structure 

of the splay soliton was calculated only by using variational approximation, (ii) The calculations were limited to 
^ , temperatures near Tc. (iii) The effect of different parameter values was not studied, (iv) The amplitudes and widths 
I of the satellite peaks as well as peaks of higher order were not studied, (v) Dissipation was neglected. It is just these 
' O points that we address in this paper. An additional motivation is that the study of solitons opens the way to detailed 
^ understanding of more complicated topological objects such as vortices. 

^ We find that, using dissipationless theory, the agreement between the theoretical and experimental frequency shifts 

• • i: is generally very good. However, we find a puzzling disagreement in the transverse oscillation frequency of the splay 
soliton. This disagreement has remained unnoticed because no detailed comparison between theory and experiment 
has been published. Further, we find that the theory is particularly inflexible to explain this discrepancy away. Taking 
^ ■ into account dissipation, in particular spin diffusion, changes these conclusions. On one hand, the disagreement in 
5^ the splay-soliton frequency is reduced. On the other hand, we find considerable extra shift of resonance frequencies 
in high fields, which has not been reported experimentally. We also point out that the longitudinal resonance of the 
splay soliton has not been studied experimentally. Measurement of these quantities would be important to test our 
understanding of the basic properties of superfluid '^He. 

We start in Sec. ^ with a short introduction to the hydrodynamic theory and NMR in ^He-A. In Sec. Illll we solve 
numerically the structure of the splay soliton. The frequency and the absorption of the principal satellite peak are 
determined in Sees. IIVI and ignoring dissipative effects. In Sec. |3 we calculate the frequencies and absorption 
of the higher modes. In Sec. IVII we take into account dissipation and calculate the effect of the spin diffusion and 
normal-superfluid relaxation on the absorption spectrum. 
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II. HYDRODYNAMIC THEORY 



Here we briefly present sorae main points of the hydrodynamic theory and NMR in "^He-A. The order parameter of 
superfluid '^He-A is a 3 x 3 tensor of the form^iii 

Af^j = Adf,{mj +inj), (1) 

where d, rh and fi are unit vectors and m _L n. It is conventional to define 1 = m x n, which gives the axis of the 
orbital angular momentum of a Cooper pair. The unit vector d defines the axis along which the spin of the Cooper 
pair vanishes. In a static magnetic field, the structure of a soliton can be determined by finding a local minimum for 
the free energySiiiii^ 

^;tatic = J dM/d+/g + /h). (2) 

Here /d comes from the dipole-dipole interaction between nuclear moments, 

/d = -iAd(d.i)2, (3) 

/h from coupling to the external field B, 

A = iAh(d.B)2, (4) 

and /g from the gradient of the order parameter, 

2/g - P±v2 + (pii - p^){l ■ v,)2 + 2Cv, • V X i 

- 2Co(i • v,)(i • V X i) + K,{V ■ i)2 + K,{1 • V X i)2 

+ i^blix (V X i)|2 + if5|(i. v)d|2 + x6^[(ix v)4f . (5) 

The gradient energy jSJ also includes the kinetic energy arising from the superfluid velocity Vg — {h/2m^) rhi^hi, 
where is the mass of a '^He atom. However, in the following we limit to the case of zero superfluid velocity. The 
parameters appearing in the gradient energy have been calculated in the weak-coupling approximation b y C rosa^ 



lation D v K^ i 

and Dorflefi^ the latter including more Fermi-liquid parameters. For numerical values see Refs. [13 andlli The 
characteristic scales are given by the dipole length = (ft/2r7i3)-^p||/Ad ~ 10 /im and the dipole field = 
^/X^ « 2 mT. 

We consider a static field Bq — Boz. We assume that the equilibrium d, denoted by do, lies in the plane perpen- 
dicular to Bq: 

do = xcos6' + ysin6'. (6) 

This situation is always achieved in large field £?o ^ i?d, where d is forced to the plane by /h Q, but in some cases 
this happens in low fields as well. Minimization of the total energy (0) gives for 9 the equation 

2?^+(i-do)(ixdo). = 0, (7) 

where the operator V is defined by 



Ad Ad 



1(1 -V)/ . (8) 



In a dynamic magnetic state one has to include the spin magnetization 7S as a new variable in addition to d and 
1. The effective energy density has the for 

/off = -S - 7S • B -)- /d -t- /g, (9) 



where 7 is the gyromagnetic ratio and X the susceptibihty tensor. This leads to the equations of motion 

S = 7SxB-dx^, (10a) 
5A 

d = d X 7 (^B - , (10b) 

where / = /d + /g and x is the susceptibihty in the normal state. The motion of 1 is strongly limited by viscosity and 
therefore we assume that 1 is independent of time>i2i Equations H10|l describe dissipationless dynamics. The inclusion 
of dissipative terms is postponed to Sec. I VII The field B is the sum of the static field Bq and a small radio- frequency 
field B' that oscillates at angular frequency lo. Throughout this paper we limit to study the linear response of S and 
d to B'. We parametrize the deviation of d with two parameters dg and dz, 

d = do + (z X do)(ie + zdz. (11) 

For S we parametrize the deviation S' from the equilibrium So — xBo/mo7 by S'^ = Sz — So and circular components 
= Sx i iSy. Similar definitions are used for other vectors as well. We linearize Eqs. pO|l and assume the time 
dependence S'(i) — S' exp(— icjt), etc. Using the equilibrium condition Q we get 

uS^ = ±uoS^TXde^''{V + Ui_)d,+iXde^^''lJ^deTxBoB±, (12a) 

^4 ^ (S-e'o - S+e-'') + 1 (B+e-'O - B.e'") , (12b) 

ujS', = -iAd(I? + C/||)de+iAd(do X i),/,4, (12c) 

2 

iude - i^S',-ijB'^. (12d) 
Here loq = jBq is the Larmor frequency. The potentials [/|| and Uj_ are defined by 

C/|l = l-ll-2i\xdo)l (13) 

= I- 2il - (i X do)^ - ^{ver . ^ef. (w) 

In order to simplify Eqs. \i2\ . let us consider the special case 1^ = 0. In this case the equations separate into 
independent blocks for longitudinal and transverse oscillations of the magnetization. The resonance frequencies are 
determined by two independent Schrodinger-type equations for de and dz, 

{V + U\\)de = a\\dg, (15) 
{V + Ui_)dz = a^dz. (16) 

The eigenvalues q;|| jj and aj^.fc of these equations are related to the resonance frequencies as 

wi^fc - ul + n^a^^k- (18) 

Here 51 = (Aio7^-^d/x)^^^ is the longitudinal resonance frequency of the A phase. The corresponding eigenfunctions of 
Eqs. l[T^ and l(TH|l are denoted by -011, fc and V'±,fc, respectively. They are assumed to be normalized: / d^r iV'H.fcP = 1 
and J d^r|i/'_L,fcP = 1- Because we have temporarily neglected dissipative processes, the power absorption P{uj) 
consists of (5-peaks, P{ijj) = X^fc Ik^iu — ojk)- 

Instead of assuming Iz = 0, an alternative approach to Eqs. 1)12(1 is to study the high-field limit ujq ^ fl. More 
precisely, one can calculate the resonance frequencies as a power series of ft^ and neglect terms of the order of 
and higher. In this approximation all the three components S"^, 5*^, and S~ decouple. The eigenvalue equations and 
frequencies are the same as above equations ()15 |l - ((18|l except that Eq. ((18|l is valid only to leading order in 51: 



fl2 ^4 -1 

^0 + 7^a±.k + 0{^) . (19) 



For the rest of this section we assume the high-field limit ujq ^ ft. 



FIG. 1: The structure of a splay soliton, where the 1 field (light arrows) has splay shape. The structure is homogeneous in the 
y-z plane of the soliton wall. Both d and 1 are perpendicular to the static field Bo = Boi. 



In the case of dipole locking, l(r) — d(r), the lowest bulk eigenvalues are a||^b = ct±.h = 1. In this case only the 
bulk eigenstate gives rise to absorption || = Vi\\ and /b,j_ = where the two modes 



(20) 
(21) 



and V = J cfir is the volume of the liquid. 

In the presence of dipole unlocking also other eigenstates contribute to the absorption. Their intensities are given 

.2.9.20 



'II, fe = «|| a||,/c Q\\,k = *l| "11, fc 



I±,k — i±Q±,k — i± 
Here the Qk's satisfy the sum rules 

oo oo 

fc=o k=a 

oo ^ 



(iVV'±,fe(r)exp[i6l(r)] 



fc=0 



'^a±^kQ±,k = 

oo „ 

7, n 



d\ [(i • do)2 - ii 

d^rUl 



fc=0 



'^oi±,kQ±,k 



k=0 



(\-dof-ll +(ixdo)2(l.do)^ 



(22) 
(23) 

(24) 
(25) 
(26) 
(27) 
(28) 



and so on. The sum rules can be derived using the orthogonality properties of the eigenfunctions. [In Eq. (|28|) one also 
needs Eq. ((T)).] The lowest-order rules apparently are equivalent to the sum rules presented by Leggett.2i 



III. EQUILIBRIUM STRUCTURE OF A SPLAY SOLITON 

The minimum of the dipole energy Q can be achieved in two ways: either 1 and d are parallel or they are antiparallel. 
This double degeneracy gives rise to domain walls or solitons. On one side of the wall 1 = d and on the other side 
1 = — d. We will study cases where the static field Bo is either perpendicular or parallel to the plane of the soliton 
wall. The soliton structures in these cases are known to have composite twist and splay structures, respectivelyi^ 
The structure and the fundamental NMR frequency of a twist soliton have been solved analyticallyiSiLS Here we 
concentrate on the case of a splay soliton, which has previously been studied only by variational methods i^i^iSSiS^ Its 
structure is illustrated in Fig. ^ 




FIG. 2: The functions ri{x) and 8{x) for a splay soliton at temperatures T/Tc — 0.01,0.25,0.50,0.75,1.00. The parameter 
values correspond to weak coupling at 3.4 MPa (Ff = 14.45 and Fi = — 1). 



In the cases of a twist or a splay soliton, both 1 and d are in the plane perpendicular to Bq. We parametrize 

1 = xcosr/ + y sin?7 (29) 

and d as above, Eq. . For a splay soliton the angles 6 and f] are taken as functions of x only. We use the boundary 
conditions ri{+oo) — 9{+oo) and ?7(— oo) = 0{—oo) + tt. Substituting the vector fields © and into the total 
energy ^ gives that the energy per unit area is 



dO 

-{Kg sin^ rj + cos^ ^ ) ^ ('^ ~ ^)] 



(30) 



The energy functional ipn|l was discretized using equal intervals {xi — —L/2 + Li/N, i — 0, . . . , N), and the values of 
r] and 6 at these discrete points were taken to be the minimizing variables. The boundary conditions were taken into 
account by linear initial approximations i] — {tt/2) — (2tt/3L)x and 9 = {tt/3L)x. The minimization can be done by 
a simple relaxation method. The resulting functions r]{x) and 0(x) at different temperatures are shown in Fig. [5] 



IV. PRINCIPAL NMR FREQUENCY 

We now apply the theory of Sec.mto calculate NMR properties of solitons. For solitons 1^ = 0, which decouples the 
longitudinal and transverse modes at arbitrary ujq in Eqs. (|12|l . Thus the eigenvalues a±^k are related to resonance 
frequencies as given in Eq. (|18|l at any field. 

The main signature of solitons in the NMR spectrum comes from the lowest eigenvalue of Schrodinger-like equations 
()15|l and 1)16(1 . This lowest frequency can be calculated, for example, using a variational formulation: 

JCelVVP + {K^ - Kg)\1 ■ VVP + AdC/|ul^p] 

" vR^W ■ ^^^^ 

This was discretized and the values of tpi^i) were taken to be the minimizing variables. The form for the initial 
approximation used was ipix) = cosh"(gx) and the boundary conditions dip/dx — were assumed at the end points. 
The length L was increased until its effect disappeared. One solution is shown in Fig. |21 In Fig. 01 we plot the 
temperature dependence of the lowest resonance frequencies as well as some experimental data. For completeness we 
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FIG. 3: The potential U± (solid line), the lowest eigenvalue a_i_, and the eigenfunction tp± (dashed line, vertical scale arbitrary) 
for transverse resonance of a splay soliton at T « Tc. The absorption thickness sj_ gives the effective width of the wave function 
in NMR. 
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FIG. 4: The reduced resonance frequencies ay 1171 1 and Ofx I18II as functions of temperature. Splay-soliton transverse resonance: 
four upper solid lines and experimental points o (Ref. ;6), A (Ref. s^nd •(> (Ref. |23). Splay-soliton longitudinal resonance: 
four lower solid lines. Twist-soliton transverse resonance: four upper dash-dotted lines and experimental points -I- (Ref. I24I1 
and □ (Ref 23). Twist-soliton longitudinal resonance: four lower dash-dotted lines and experimental points x (Ref. |a) and V 
(Ref. 0). In each case there are four theoretical lines, which in order of decreasing a correspond to = & weak coupling, 
-F" = & gap enhancement, = — 1 & weak coupling, and F^ = — 1 & gap enhancement. Good agreement between theory 
and experiment is achieved with F^ = — 1 and weak coupling except for the high-temperature transverse splay data. The effect 
of strong-coupling corrections is shown at T/Tc = 1 where the lower ends of the bars correspond to 7 = 2 and the upper ends 
to 7 = 4. The experimental data is plotted with higher resolution in Fig. 1101 



also include the analytical results for the twist solitontSi^ 



twist _ -^6 

The reduced frequency shifts a\\ H17(l and a± (|18(l depend only on the ratios of the hydrodynamic coefficients K^, 
Ks, Kt, K5, and Kq. (Note that the absolute magnitudes of Ki's and Ad define length and energy scales that do not 
affect q;|| or aj_.) In the weak-coupling approximation the ratios of if^'s are functions of an infinite set of Fermi-liquid 



(32) 
(33) 



parameters and F"-, with j = 1, 3, 5, etcii^ Here we neglect all the coefScients with j > 1 since they are unknown. 
For F( we use the value by Greywall^^ at the melting pressure, Ff — 14.5. However, the dependence on Ff is weak. 
For example, the variation Ff = 14.5 ± 1 shows up only in the fourth decimal of a at T = 0.8Tc. In terms of pressure, 
the maximum difference in a between 2.6 and 3.4 MPa is 1% at temperatures higher than 0.5Tc. The value = —1 
was taken from Ref. and it is also consistent with Ref. In order to see the effect of -Ff we also used F^ 0. It 
can be seen in Fig. ^ that this shifts the resonance frequency up at temperatures below T^. 

There are no quantitative calculations of strong-coupling effects in the A phase at a general temperature. In order 
to get some idea how strong coupling could affect the soliton frequencies, we use a "trivial strong-coupling" model 
developed by Serene and Rainer^'* for the B phase. In this model the weak-coupling energy gap is multiplied by a 
factor that depends on the temperature and on the jump of the specific heat ACe/Cn at T = Tc. This dependence 
is tabulated in Ref. We adapt this model to the A phase by calculating the multiplying factor using the same 
table but substituting |ACA/Cn in place of ACe/Cn, and using extrapolation when needed. We take ACa/Cii from 
measurements by Greywall.^^ It can be seen in Fig. ^ that the gap enhancement, which affects only intermediate 
temperatures, has a smaller effect than the change of Ff . 

In the limiting case T Tc the reduced frequencies a|| and a± are independent of any parameters appearing in the 
weak-coupling model, including also the gap enhancement. This is a consequence of the Ginzburg-Landau expansion 
that gives to the parameters Ki the ratios i^b : Kg : Kt : : Kq = 7:1:1:2: (7 -f-l). The only free parameter here 
is 7, which in the weak coupling (with or without gap enhancement) has the value 7 = 3. This value is changed only 
when nontrivial strong-coupling corrections are included. The "weak-coupling-plus" model by Serene and Raines^S 
gives an estimate 7 ~ 3.12. Figure 01 shows the reduced frequencies corresponding to 7 = 2 and 7 = 4. We see that 
this variation changes the splay soliton a± less than 3% but for the twist soliton the effect is ten times larger. 

There is rather good agreement between the experiment and the theory corresponding to weak coupling and 
Ff = —1. Equally good agreement is achieved with gap enhancement and « —0.7. The longitudinal twist data 
is a strong indication that the deviation from the weak-coupling value 7 = 3 is small, as predicted by Serene and 
RainerS The only major difference between theory and experiment exists in the transverse splay-soliton frequency 
at high temperatures. It seems very difficult to improve the agreement in this case by any change of the parameters 
in the theory above since a^''^'*' at is effectively fixed. On the experimental side, one possibility is that the field is 
not precisely in the plane of the soliton in the measurements. This would add a small twist component to the splay 
soliton and thus shift up the frequencyi^ Another possibility is that the inclusion of relaxation mechanisms could shift 
the calculated resonance frequency, as will be discussed in Sec. I VII 

At Tc we find the eigenvalues ay = 0.388 and a± = 0.672 for 7 = 3. These differ slightly from the variational 
results by Maki and Kumar^ that are 0.403 and 0.677, respectively. 



V. ABSORPTION AND HIGHER MODES 



Here we calculate the intensity of the principal soliton peak and analyze the absorption at other frequencies. For 
planar objects intensities (|22|l and (|23|l are most conveniently expressed in the form of an absorption thickness Sk, 
which equals Qk divided by the area of the planar object, 

= ^, (34) 

for each mode k. The absorption thicknesses in the lowest eigenstate are plotted in Fig. [S] 

The absorption for a twist soliton can be calculated analytically. Substitutingthe eigenfunction tp = ci sech^ az 
into Eqs. I|22|l and 1)23(1 and using 9 — C2 ~ sgn(z) arccos(sech az) (Refs. and^ we find 

II „r(^||)r2(^) 
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FIG. 5: Absorption l,'i4^ of the principal soliton peak in four cases. The parameter values correspond to weak coupUng at 
3.4 MPa {F( = 14.45 and = -1). The numerical values at Tc are sj^''^ = 6.63Cd, s'f''^ = 7.35Cd, Sy"''"* = 8.64^d, and 



Here T is the Gamma function and 

CL = \ 

1 / j mt + K^ \ 

In order to get all the resonance modes, we discretize the space as discussed above [xi = — L/2 + Li/7V, i = 0, . . . , iV, 
for twist soliton replace x hy z). This means that Eqs. p5|) and turn into a matrix eigenvalue problem. This 
can be solved by standard library routines for matrices of reasonable size. This gives all the frequencies, including the 
fundamental one However, the higher modes extend over the whole interval, and it is important to fix boundary 
conditions for them. We require zero derivative oi ij] &t x — ±L/2. This can be justified by considering a lattice of 
solitons. The general eigcnfunctions of Eqs. \ib\ and are of the Bloch form, but only strictly periodic functions 
lead to nonzero absorption in Eqs. H22I) and H23() . The unit cell of a soliton lattice x = — L/2, . . . , 3L/2 consists of 
two solitons located around a; = and x — L. The vectors 1 and d have the symmetries \{L/2 + x) — \{L/2 — x) 
and d(_L/2 + x) = d(L/2 — x). This implies that the eigcnfunctions can be classified as symmetric or antisymmetric 
(with respect to a; = L/2). The symmetric solutions have dip/dx{±L/2) — 0, and the antisymmetric solutions can be 
neglected since they do not contribute to the absorption. 

The results for the frequencies and intensities of the higher modes are shown in Fig. As expected, the higher 
modes depend on the size L of the system used in the calculation. The modes seem to appear as pairs above the 
bulk peak. The bulk peak also seems to consist of two peaks that do not appear exactly at a = 1 due to finite L. In 
fact, the finite L corresponds to a periodic lattice of solitons, which should give rise to peaks that are analogous to 
the Bragg peaks in the x-ray scattering from periodic solids. Such peaks are indeed seen experimentally in the case 
of a vortex sheet, where a nearly periodic arrangement is automatically generated^ The observation of the peaks 
depends essentially on dissipative effects, which broadens the peaks. This will be discussed in the following section. 



(37) 

(38) 
(39) 



VI. THE EFFECT OF DISSIPATION 



There are two sources of dissipation that are important for solitons. The nonequilibrium between normal and 
superfiuid magnetizations is the main relaxation mechanism in homogeneous superfluid ^jjg 3i,3,2.,a3_ j-,^ inhomogeneous 
situations the nonequilibrium between normal magnetizations at different locations causes relaxation via spin diffusion. 
We study these two mechanisms below, the latter in the case where the inhomogeneity arises from the soliton. We 
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FIG. 6; The higher modes in the case of transverse resonance of a splay soliton. The length used in the calculation is L = 50^d 
in the main frame whereas the inset also contains the points for L = 250^d- In both cases the lattice spacing of the discretized 
lattice is O.l^d. Other parameters are weak coupling at 3.4 MPa (Ff = 14.45, = — 1) and temperature T = 0.95Tc. 



neglect extrinsic effects such as the inhomogeneity caused by nonuniform magnetic field. We also do not consider 
dissipation caused by orbital motion. According to the estimation in Ref. Ilfll this is important only at temperatures 
very near Tc, region which is vanishingly small on the scale of Fig.^ 

Our treatment of the dissipation is purely phenomenological. For normal-superfluid relaxation we use the model 
developed by Leggett and Takagi, which is extensively described in Ref. 33. For the relaxation of spatially noncqui- 
librium magnetization we use simple diffusion equation. This approach can be correct only when the mean free path 
of the quasiparticles is smaller than the thickness of the soliton, and therefore it necessarily fails at low temperatures, 
where the mean free path diverges. Our treatment also neglects the tensor character of the spin-diffusion coefficient, 
as well as the anisotropic energy gap with varying anisotropy axis l(r), which leads to Andreev reflection of the 
quasiparticles. 

In order to treat the relaxation, the spin polarization S is divided into a superfluid (pair) part Sp and a normal 
(quasiparticle) part Sg = S — Sp. Both parts have separate equations of motion 



7S, X (^B - /^o7^Sp) 



7Sp X B 



MoT — 



7d X 



B 



Mo7- 



Xo ' 



Xo 



A)Sp 

A)Sp 
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Axo 



AS,] + kV^S, 



AS„ 



d X — 
dd 



(40a) 
(40b) 
(40c) 



Except the spin-diffusion term kV^S,, these equations are the same as in Ref. 's^ Here xo is the susceptibility in the 
absence of Fermi-liquid effects. As above, x is the susceptibility in the normal state, so that x/xo = 1/(1 + Fo)- "^^^ 
function X{T/Tc) is defined as the equilibrium fraction of the superfluid magnetization, and is given by Eq. (4.23) in 
Ref.'s^. The Leggett- Takagi relaxation time r describes local relaxation of Sp and S, to their equilibrium values. The 
spin diffusion appears via term kV^S, in the equation for Sq (|40a|l . In the normal state the spin-diffusion constant 
K is related to the spin-diffusion time td by^ 



(41) 



where vf is the Fermi velocity. We use this parametrization also in the superfluid state. Note that k is defined as the 
spin-diffusion constant for Sq, so that the effective diffusion constant for the total magnetization S is (1 — X)k. 
We continue by writing the equations for variables d, S, and the deviation from local equilibrium 



r? = Sp - AS. 



(42) 



Similarly, as in Sec.^ we linearize the equations and assume harmonic time dependence. We use the fact that Iz — 
in solitons. As a consequence, the transverse and longitudinal modes separate at any field Bq. Differing from Sec. 



introduce dimensionless quantities by defining S"^ — {n"^ / XdOJo){S'^ ± iSy) and 77=*= — {fl'^ / XdiLJo){Vx ^ i^y)- For 
transverse oscillations we get the equations 
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The full solution of this problem can be written as a five-component vector ip = {S^ rj^ dz tj^ S*^)^ = '^^Ciipi, 
where the i/^i's satisfy the homogeneous equation where B± — 0. The eigenvalues w^ of the homogeneous equation are 
now complex valued, and the eigenvectors ipi are not orthogonal to each other. In order to solve for the absorption 
spectrum we find that the following adjoint eigenvalue problem:^^ 
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(44a) 
(44b) 
(44c) 



has the same eigenvalues w^ as the homogeneous version of problem (|43|l and that the eigenvectors V'j = 
{S" rj^ d^A 1]^ are orthogonal to ipj's: 



(45) 



When deriving this we have assumed zero derivate for the wave functions far from the soliton by considering a lattice 
of solitons as was discussed in Sec. The power absorption can now be written as 



P(w) = ]puj j Im(B'-S')dV 
^Im j 



4^*0 



d:^r{B*_S- 



(46) 



where the coefficients Cj are given by 
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B_S- - B+Sj + -(e-'^B+ - e'^B_)d^ 



(47) 



We observe that the eigenfunctions are not symmetric with respect to the center of the soliton. This apparently is 
caused by the spin-diffusion term. 

In the hydrodynamic limit where wr ^ 1 one may solve 77^ from Eq. H43|l and to linear order in r 



ri^ = -tA[L - AjKV"^^ ± iT[L - A) — e^-"[V + U )dz. (48) 

Wo 

If we additionally ignore the term \/^ri^, the equations for the three-component vector {S^ d^ S^)'^ read 

K 



w 
Wo 
w 



±1 + i(l - A) — V 

Wo 



s- 



Wo 



Xo 

.Txn^ii-\),^ .... 1 



1 



Bo' 
Xo 



s+ 



iV + U^)dz + — {B+e-''-B^e'') 
Xo^oA 2Bo 



(49a) 



(49b) 



The homogeneous adjoint problem for this is given by 



Wo 



±l + i(l-A) — 



XO 



Wo 



UJn 



Xo^oA 



(50a) 
(50b) 



The equation for the power absorption does not need any modifications. One must be careful since the validity region 
of the hydrodynamic approximation is not very large. For temperatures near Tc one must, for typical values of t, have 
-Bo ^ 15 mT. We mainly used the hydrodynamic approximation for checking our calculations in the limit of t ^ 0. 

For longitudinal case, where B' — B'z \\ Bo, we also define dimcnsionless quantities by writing S' = zAd5z/fi, 
ij' — zAar/z/f^, and d' — dgz x do. The equations of motions are 



^S, = i(l-A)-V25,-i-V277,-i(I? + C/||)de, 



= -iA(l-A)^V^^.-i(-L 



A-V2 ) ,7,-1(1 -A) (I? + C/||)de 



and the homogeneous adjoint problem for this is given by 



^i(l-A)iiv^&-iA(l-A)iiv',^ 



1 ^ «^^9 



It — dg 



i(I? + C/||) ^, + (l-A)ry^ 



Power absorption for longitudinal case is given by 



P(w) 



(^J / 



(51a) 
(51b) 

(51c) 

(52a) 
(52b) 
(52c) 

(53) 



where the c^'s are the coefficients of the full solution -0 = {Sz tJz dg)'^ = ^ Cj{Sz.j tJzj dgj)'^ and given by 



da id^r 



(54) 



Similar to the transverse case one could write the hydrodynamic equations using only dg and Sz ■ 

The numerical solution for these different eigenvalue problems is obtained by dividing the calculation length x ~ 
~L/2, . . . , L/2, for example, to 1000 points and approximating the spatial derivatives by differences. The eigenvalues 
and vectors of the resulting sparse 5000 x 5000 (or in the longitudinal case 3000 x 3000) matrix A are then solved by 
MATLAB. We make use of sparse matrices and calculate normally only 20 lowest eigenvalues with Re[a;fe] > that 
give the main contribution to the absorption spectrum. The complex mode frequencies LUk are related to the reduced 
frequency shifts ak as given by Eqs. (fTTjl and (fH^ . 

Additional parameters appear in the calculation compared to the dissipationless case. The most crucial ones are the 
relaxation times r and td . The Leggett-Takagi time t can be extracted from the width of the bulk peak. Measurements 
of longitudinal and transverse resonances are in good agreementf^S, We have reanalyzed the data by Gully et al.^ 
including strong-coupling corrections as described in Sec. II VI and get the fit r = [1.70 + 6.71(1 — T/Tc)]{l + Fg )10^'' 
s in the range of T/Tc — 0.78-0.98. We fi x Fn = —0.746. For the spin-diffusion coefficient k we use the value 
kT^ = 1.1 X 10"^ mK^/s from Refs. [slandU^ This corresponds to the spin-diffusion time td = 1.25 x 10~''T~2 
mK^ s. A few absorption spectra are plotted in Figs. CHHI 

The effect of the normal-superfluid relaxation seems to be simple line broadening. In the transverse resonance the 
linewidth measured on the a± scale is inversely proportional to the magnetic field Bq for Bq ^ B^. This is caused by 
i/(a;oT) term in Eq. (|43b|l . There is no field dependence in the longitudinal resonance since Eqs. t\51\\ do not contain 
loq. In the longitudinal case the linewidth on the a scale has no strong temperature dependence and approaches a 
finite constant when T ^ Tc- In the transverse case the linewidth on the a scale vanishes as T — > Tc- AH these 
characteristics are the same as predicted for the bulk peaks. 




FIG. 7: Absorption spectra for transverse twist soliton (dash-dotted lines) and splay soliton (solid lines) at T/Tc = 0.92, 0.95, 
and 0.97 when L = 50Cd, Bo = 3 mT, r = [1.70 + 6.71(1 - T /Tc)\(l + Fo'')10"^ s, and td = 1.25 x lO'^'T"^ mK^ s. Other 
parameters are obtained using weak coupling and pressure of 3.4 MPa with = —1 and Fq — —0.746. 




FIG. 8: Transverse absorption spectrum for a twist soliton at T = 0.75Tc and Bo ~ 28.4 mT (solid line). Other parameters are 
the same as in Fig. |7| The effect of the two dissipation mechanisms is demonstrated by setting r = (dash-dotted line) and 
Td = (dashed line). 



The effect of spin diffusion is more complicated. In addition to broadening, it shifts the principal soliton peak 
to higher frequencies. Since the spin diffusion is due to inhomogeneous texture, its strength in the transverse case 
is obtained by comparing the term i(l — X){k/ujo)'^^S^ with the term {Q'^ /ujq) exp{±\6){T> + U±)dz in Eq. (|4;^a|l . 
Therefore, the effective diffusion constant K^g = (1 — X)KU!o/fl^ is linearly proportional to the magnetic field. As a 
consequence also the reduced frequency a± of the satellite peak is field dependent. When Bq = 5 mT the effect of 
diffusion is already quite large at T = 0.95Tc and the soliton satellite peak is almost smeared out. In the longitudinal 
case there is no field dependence and the effective diffusion constant reduces to K^g = {1 — X)K/rt. In both transverse 

and longitudinal cases the effective diffusion coefficients K^^g and K^g diverge when T ^ T^, as 17 ^ 0. 

The relative contributions of normal-superfluid relaxation and spin diffusion are illustrated in Figs. |H1 and El where 
both are separately shut off by setting r = or td =0. It can be seen that the broadening of the longitudinal satellite 
peak is mostly caused by normal-superfluid conversion at T = O.QSTc- In the transverse case at high field the spin 
diffusion becomes the dominant relaxation mechanism, as well as approaching Tc. In uniform order parameter the 
width of the bulk peak comes solely from normal-superfluid relaxation. Our soliton lattice has many higher peaks in 




FIG. 9: Longitudinal absorption spectrum for a splay soliton at T = 0.95Tc (solid line). There is no dependence on the magnetic 
field, and other parameters are the same as in Fig.|3 The eflPect of the two dissipation mechanisms is demonstrated by setting 
r = (dash-dotted line) and to ~ (dashed line). 




FIG. 10: The resonance frequencies and linewidths as a function of temperature at Bo = 3 mT. The solid lines give the frequency 
of maximum satellite absorption, the dash-dotted lines 75%, and dotted lines 50% of the maximum. The dashed lines give 
resonance frequency in the absence of dissipation. Other parameters are L = 50^d, t = [1.70 -I- 6.71(1 — T/Tc)]{l + Fo)W~^ s, 
TO = 1.25 X 10"^r-2 mK^ s, pressure 3.4 MPa (Ff" = 14.45, = -1, and Fg = -0.746 and weak coupling). The experimental 
points are the same as in Fig. |1] 



the absence of dissipation, but dissipation seems to remove them. Note that the higher peaks are mostly suppressed 
by spin diffusion in the case of Fig. |S1 but by normal-superfluid relaxation in the case of Fig. O 

The relatively high field in Fig. [S] has been chosen to allow a comparison to the spectrum of Hakonen et alm^ The 
shape of the twist-soliton spectra are very similar even if the temperature is quite low, where our modeling of the spin 
diffusion is under suspect. Further spectra at lower field are shown by Parts et alw^ The comparison to Ref. ^is a 
different case. In this experiment such a dense soliton lattice was created that no bulk peak was observable. Avenel 
et al^ and Gould and Lee^ show soliton spectra measured by sweeping temperature. In all cases a proper comparison 
to our calculations would require first estimation of the density of solitons and then correcting our calculations for 
that density, which we have not done. What we can state that there seems to be no obvious contradictions between 
theory and experiment concerning the linewidth. 

A summary of our results at i?o = 3 mT is shown Fig. 1101 We see that the frequency of maximum absorption is 
shifted towards higher frequencies, especially at temperatures near T^. This alleviates the disagreement found between 
dissipationless theory and the experiment by Gould and Lee^ in the transverse mode of splay soliton. However, the 



experimental data does not show any sign of the divergence at predicted by our modeL The shift of the resonance 
frequency results from spin diffusion. For the longitudinal mode it is field independent, but for the transverse mode 
it increases with increasing field. In high field this shift can be substantial even at T « 0.7Tc, see Fig. |H1 This 
poses a problem since no field-dependent shift has been observed experimentally. This is the case particularly for the 
transverse mode of twist soliton measured at Bq = 9.9 mT (□) (Ref. l23j) and 28 mT (+) (Ref. l24j) . Similar evidence 
for transverse mode of splay soliton measured at 2.5 . . . 3.7 mT (o) (Ref. 0) and at 15 ... 20 mT (A) (Ref. Ii(T[i is not 
so clear since the latter data is measured under different conditions, as discussed above. 



VII. CONCLUSIONS 



The frequencies of soliton satellite peaks in the NMR spectrum are calculated at all temperatures. The agreement 
of dissipationless theory with experiments is very good. However, there is a small difference in the transverse-mode 
frequency of a splay soliton near Tc. This difference is partly explained by taking into account spin diffusion. The 
spin diffusion shifts up the reduced frequencies a at high fields, which has not been observed experimentally. We 
hope that new experiments could clarify this problem. We also point out the relatively narrow longitudinal line of 
the splay soliton, which has not been studied experimentally. 

It would be of interest to extend the present calculations to vortex lines in '^He-A. This requires two-dimensional 
calculation and thus would be computationally much more demanding than the present one. 
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